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SUMMARY 

The general-purpose finite-element program MSC/NASTRAN Is used to study 
the Interply layer effects on the free-edge stress field of symmetric angle- 
plied laminates subjected to uniform tensile stress. The free-edge region Is 
modeled as a separate substructure (superelement) which enables easy mesh re- 
finement and provides the flexibility to move the superelement along the edge. 
The results Indicate that the Interply layer reduces the stress Intensity sig- 
nificantly at the free edge. Another Important observation of the study Is 
that the failures observed near free edges of these types of laminates could 
have been caused by the Interlaminar shear stresses. 


INTRODUCTION 
Preliminary Remarks 

Free edges In composite laminates often promote complex three-dimensional 
stress fields with steep gradients. This phenomenon has been studied both 
experimentally and analytically to a considerable depth during the last 15 
years. It Is thought that delamlnatlons and matrix microcracking occurring 
near free edges are due to the steep Interlaminar stress gradients In the free- 
edge region. However, much controversy still exists regarding the nature and 
magnitude of the Interlaminar free-edge stresses. The free-edge region of a 
composite laminate Is usually modeled as homogeneous anisotropic layers sharing 
common boundaries between them. In reality, however, a thin, soft matrix layer 
exists between the plies. The effect of such a layer on the free-edge stress 
field must be considered In order to completely describe the stress field as 
the free edge Is approached and also to better understand the apparent singular 
behavior of some of the stresses near a free edge. The singularities are In- 
herent to the formulations when there Is an abrupt change of material proper- 
ties across the common boundary as when the Interply layer Is neglected. 
Furthermore, singularity Is only a mathematical concept and not a realistic 
Issue In practical problems. 


Background 

Numerous papers pertaining to the free-edge stress have been published 
during the last 15 years. A brief description of the literature on the subject 
Is given In order to provide the background and Impetus leading to this 
Investigation. 


*NRC-NASA Research Associate. 


The problem of finite-width laminate under uniform extenslonal strain was 
first analyzed by Pipes and Pagano (ref. 1). The formulation was based on 
anisotropic elasticity, and the solution was obtained using a finite-difference 
scheme. A number of Investigations have subsequently been conducted on this 
problem with a fairly wide range of approaches. Quasl-three-dlmenslonal 
finite-element approaches which use plane-strain plate elements can be found 
Iri references 2 to 4. Hsu and Herakovich (ref. 6) obtained a solution using a 
perturbation technique. Singular hybrid-stress finite-element models were 
employed by Spllker and Chou (ref. 6) and Wang and Yuan (ref. 7). The latter 
work Is based on the theory of anisotropic elasticity and Lekhnl tskl 1 1 s stress 
potentials (ref. 8). A fully three-dimensional formulation with a three- 
dimensional finite-difference solution scheme Is given In reference 9. Re- 
cently, a global-local laminate variational model has been proposed by Sonl 
and Pagano (ref. 10). A review of the free-edge stress problem with a compar- 
ative assessment can be found In references 4 and 11 to 13. 

Approximate procedures, which assume specific form for the stress dis- 
tributions near the free edge and use elastic equations of equilibrium, have 
also been employed with some success (refs. 14 to 16). The main Intention of 
these works Is to present procedures to design laminate with appropriate 
stacking sequence such that the resulting laminate Is less likely to develop 
high Interlaminar stresses and subsequent failures In the free-edge region. 


Synopsis 

A survey of the literature on the free-edge stress problems leads to the 
following Important concluslons/lssues: 

(1) A fully, three-dimensional finite-element analysis or elasticity 
solution Is not available at the present time. 

(2) The Interply layer of finite thickness has not been considered In 
the solution available to date. 

(3) -Much controversy exists regarding the nature and magnitude of the 
Interlaminar stresses near the free edge. 

(4) The free-edge region has been customarily modeled as dissimilar an- 
isotropic layers sharing common Interfaces. Consequently, stress singularities 
are Inherent to such formulations. 


Present Investigation 

This investigation Is directed toward addressing the previous Issues with 
the aid of a finite-element analysis. In a typical solution scheme of finite- 
element or finite-difference type, the layers are divided into a fine mesh. 

The mesh size near the Interface Is usually of the order of the dimensions of 
the fiber and the Interfiber spacing. Therefore, a more realistic approach 
from the micromechanics viewpoint would be to provide a thin Interply layer at 
the Interface as Illustrated In figure 1. In this analysis, a thin, soft. 
Isotropic matrix layer Is provided at the Interface. This layer is referred 
to as the Interply layer In this discussion. A fully, three-dimensional 
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finite-element analysis Is carried out using 20-node Isoparametric brick ele- 
ments that are available In the general-purpose finite-element code 
MSC/NASTRAN. Another feature unique to the present study Is the use of a spe- 
cial free-edge superelement which allows the free-edge region to be analyzed 
as a substructure. This feature has the added advantage of remeshing and re- 
analyzing the free-edge region to the extent of detail desired without solving 
the primary structure every time. The analysis Is conducted on the complete 
structure. The dimensions of the laminate are chosen to reflect a realistic 
plate type structure (28 In. by 7 In. by 0.0202 In.). 


ANALYSIS 

In this section, the specimen geometry, the finite-element Idealization 
of the primary and the superelement structures, and the cases Investigated are 
described. 


Specimen Geometry 

The geometry of the laminated plate under study Is shown In figure 2. As 
shown In the figure, the length of the plate Is 28 In., the width 7 In., and 
the thickness 0.0202 in. The plate Is composed of seven layers (four plies 
and three Interply layers). The plies are arranged symmetrically to form a 
(±e) s type angleplled laminate. The coordinate axes are as shown. 


Finite-Element Modeling 

Figures 3 and A give the details of the finite-element mesh employed In 
the present analysis for the primary and the superelement structures, respec- 
tively. One element Is used through the thickness of each layer. There are 7 
elements across the plate width and 28 elements along the plate length. A 
total of 1365 brick elements are used to model the primary structure. The 
superelement contains 224 elements with progressively decreasing mesh size 
toward the free edge. With this type of arrangement, the center of an element 
nearest to the edge Is 1/512 In. from the edge. The boundary conditions Im- 
posed are such that along the edge x = 0 every node in the midplane Is con- 
strained from displacement In x-, y-, and z-dlrectlons . Also, the rest of 
the nodes along the edge x = 0 are constrained from displacement In the 
x-dlrectlon. Along the edge x = 28 in. a uniform stress of magnitude unity 
Is prescribed. 


Cases Investigated 

The types of laminates Investigated in the present study are (±10) s 
through (±80) s at 10° intervals. The (±45) s laminate Is also studied as a 
special case since It has been given substantial attention In the literature. 
The Interply layers are assumed to be Isotropic with E * 0.5 mpsl and 
v = 0.35. These are representative properties of an epoxy matrix. The thick- 
ness of the interply layer Is taken as 0.000067 In. The plies are assumed to 
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be homogeneous and orthotropic. The material properties chosen for the plies 
are typical of an AS-graphlte fiber/epoxy composite system. The stress-strain 
relation for the material Is given by 
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The NASTRAN Input data material properties for a 0°, 0.005-1n. -thick ply 
are as follows (units are In mps 1 ) : 


Gn b 19.4 

G12 = 0,5164 

G 13 = 0.5164 

g 22 

= 1.384 

G 23 = 0.4473 

G33 = 1.483 

G 44 = 0.6667 

g 55 

= 0.3974 


655 = 0.6667 

The rest of the G values are zero. 

These properties are calculated using ICAN, a resident composite mechanics 
computer program with a data base of fiber and matrix properties (ref. 17). In 
addition 'to the previous cases, a limited study Is also conducted where the 
thickness of the plate Is varied from 1/7 through 1/40 of the width. 


RESULTS AND DISCUSSION 
Free-Edge Stresses 

In order to study the three-dimensional stress behavior near the straight 
free edge, a small region close to edge Is Isolated and modeled as a substruc- 
ture (superelement). In this region, the element size Is decreased progres- 
sively as the free edge Is approached. With this type of mesh, Increasing 
levels of mesh refinement can be obtained as Illustrated In figure 5. The 
results shown In this section were obtained by employing a mesh refinement 
such that the center of the element closest to the free edge was 1/512 In. 
from the free edge. The six stresses, which were obtained near the free-edge 
region, are summarized In figures 6 to 14. These results are plotted against 
a length parameter Xf along a central line perpendicular to the free edge. 

The length parameter Xf Is related to the distance from free edge Lf by 

X f - 1 - L f (2) 
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In these figure the dashed lines represent the stress computed at the center 
of (+0) ply elements. The solid lines represent the stress computed at the 
center of the Interply layer between (+0) and (- 0 ) layers. The stresses 
are normalized with respect to the applied stress cr XXoo . The classical 
laminate analysis predictions for these normalized stresses are 


xx /c, xx« 

( 3 ) 

/ o - 

yy xx» 


The Inplane shear stress cr X y Is a function of ply angle 0 . The vari- 
ation of o X y In the to ply as a function of ply orientation Is shown In 
figure 15. The general behavior of the stress field as the free edge Is 
approached Is now summarized. 


Inplane stresses o x y . and o »y. - The Inplane stresses cr xx , cr X y, 
and oyy approach classical laminate theory predictions In the Interior of 
the ply. Close to the free edge, however, significant changes In these 
stresses are clearly seen. The maximum change In the magnitude of <j xx Is 
30 percent greater than that predicted by theory (CLT). 

The transverse normal stress oyy In the +0 layer appears to be 
negligible. However, It does not vanish because of the presence of the Inter- 
ply layer. The maximum «jyy observed In the +0 layer Is less than 10 
percent. In the Interply layer, <jy y Is significant for layups (±30) s to 
(±80) s with magnitudes ranging from 10 to 20 percent of the applied stress. 
Also, the sign of <?yy changes for (±70) s and (±80) s laminates. 

The c X y stress In the Interply layer Is negligible and can not be 
seen distinctly In figures 6 to 14. Both a xy and d xx In the +0 
layer are seen to approach a finite value near the edge. A sign reversal Is 
observed for these stresses In the case of ( ±1 0 ) s laminate. 


Interlaminar stresses j y ; . a l u and <g 7 ». ~ The Interlaminar shear 
stress a xz has the greatest magnitude among these stresses In the laminates 
range (±30) s through (±60) s . A steep Increase In magnitude Is seen near 
the free-edge zone. The stress magnitude In the Interply Is higher than that 
of the +0 layer. A maximum of 28 percent of the applied stress Is seen In 
the Interply layer of the (±40) s laminate. 

The Interlaminar normal stress a zz appears to be of minor significance. 
However, It does have a definite trend. Two sign reversals are clearly noted. 
It starts as compressive stress In the Interior of the ply or the Interply 
layer, becomes tensile stress for a short distance closer to the edge, and then 
reverts to compressive stress near the edge. The maximum <j zz observed Is 
about 3 percent of the applied stress. 

On the basis of these results, H can be concluded that <j xz Is the 
only Interlaminar stress that could Initiate delamlnatlon at a free edge In 
angleplled laminates, specifically In the sensitive range of 0 = 30° through 
0 = 60°. The Interlaminar stress <? zz by virtue of Its magnitude and 
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sign (compressive In a major portion of the free edge region) does not appear 
likely to Initiate peeling-off-type delamlnatlons near the free edge. The 
Interlaminar shear stress o 2 y Is the least significant of these stresses 
throughout the range considered. It can be neglected for all practical 
purposes. 


Effect of Ply Angle on Free-Edge Stresses 

The peak values (computed 1/512 In. from the free edge) of the Interlami- 
nar stresses a xz and a zz and the inplane stress <jyy In the 1 nterply 
layer between +0 and -e plies are shown graphically in figures 16 to 18 
plotted against the laminate type. As seen In figure 16, a xz appears to be 
very sensitive in the range (+30) s through (£60) s . The trend appears to be 
similar to the functional behavior of <* X y stress shown In figure 15. The 
maximum a xz Is observed for the (+45) s laminate (about 28 percent of the 
applied stress). 

The «j 2Z peak distribution Indicates that for a large range of e's 
the stress Is compressive. The dashed line In the figure refers to stress 
computations 3/512 In. from the free edge. The solid line refers to stress 
computations 1/512 In. from the free edge. The maximum magnitude of or 2Z 
(about 3 percent of the applied stress) Is observed for (±70) s laminate. 

The cfyy stress appears to have a definite trend with maxima at (±40) s 
and (±80) s . The maximum value of oyy Is about 15 percent of the applied 
stress. 


Influence of Progressive Mesh Refinement 
on Interlaminar Stress Distribution 

A (±30) symmetric laminate Is chosen for this study. The mesh progressive 
refinement shown in figure 5 Is employed to achieve three levels of fineness. 
The centers of the element closest to the free edge are 1/32, 1/64, and 1/128 
In. from the edge for levels 1, 2, and 3, respectively. The Interlaminar 
stresses a zz and a Z y are selected for the following discussion, and the 
results appear In figures 19 and 20. Both stress distributions Indicate sta- 
ble, converged behavior except for the region close to the edge. The Inter- 
laminar shear stress a™ distribution, however, shows much higher noise 
near the edge though the magnitudes are negligible. It Is known from the lam- 
inate equilibrium consideration that <j Z y should be zero on the edge. 


Effect of Width to Thickness Ratio on Interlaminar Stress Peaks 

A limited study with (±45) symmetric laminate has been conducted to assess 
the effect of width to thickness ratio on the Interlaminar stress peaks. The 
stresses a xz and a zz are selected for the study because they appeared 
to be the only significant Interlaminar stresses. The results are shown In 
figure 21. The stresses are computed 3/512 In. from the free edge. Four ra- 
tios of W/h, where W Is the width and h the thickness, are chosen In the 
study. The results show significant sensitivity to the W/h ratio. Both 
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stresses appear to have substantial magnitudes for the laminate with W/h * 20. 
As the laminate becomes thinner, however, o 22 appears to reach a zero value 
while o X2 continues to have a significant magnitude. 


SUMMARY OF RESULTS 

The significant results from the present Investigation are now summarized, 

1. The major stresses In the free-edge region are the Inplane stresses 
o xx , Oyy , and <# xy . The Inplane stresses cr xx and o xy Increase 

above the CLT predicted values by as much as 30 percent of the applied stress 
In the +e layers. The Inplane stress o y y Is nontrivial because of the 
presence of the Interply layer. A maximum of 15- percent of the applied stress 
Is observed In the layers of (±40) s and (+80) s laminates. 

2. The Interlaminar stresses <j x2 , a zz , and a z y are of second order 
compared to o xx In general. However, the Interlaminar shear stress <j X 2 
could become substantial In the Interply layer free-edge region depending on 
the configuration of the angleply laminate. The range of laminates from 
(±30)c through ( ±60) s are likely to develop relatively high o X2 stress 
magnitudes that may Initiate edge delamlnatlons. 

3. The peeling off stress cr 22 Is not significant, and It Is not 
likely to Initiate edge delamlnatlons In ( ±0) s angleplled laminates. 

4. Progressive mesh refinement is an effective way to describe the be- 
havior of the Interlaminar stresses as the free edge Is approached. The re- 
sults Indicate stable convergence of these stresses with some noise as the 
free edge Is approached. 

5. The magnitudes of the Interlaminar stress peaks are sensitive to the 
width to thickness ratio of the laminate. Results obtained from relatively 
narrow specimens are not indicative of practical laminate behavior. 

6. The Interply layer appears to alter the apparent singular behavior of 
the free-edge stresses. The present study Indicates finite magnitudes for all 
the stresses as the free edge Is approached. Another effect of the Interply 
layer Is to Induce transverse normal stresses cr yy In the laminate. 
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Figure 2. - Symmetric angleplied laminate of (±e) type with Interply layers 
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Figure 6. - Free-edge stresses in +10° and interply layers of (±10) s laminate. (Solid curve 
interply layer; dashed curve, center of ply.) 
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Figure 12. - Free-edge stresses in +60° and interply layers of (* 60)3 laminate. (Solid curve, 
interply layer; dashed curve, center of ply.) 
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Figure 15. - Variation of inplane stresses with ply angle 
predicted by classical laminate theory. 















Figure 21 . - Influence of width to thickness ratio on interlaminar normal and 
shear stress behavior at free edge in interply layer of (*45)5 laminate. 


